Mapping donor electron wavefunction deformations at sub-Bohr orbit resolution 
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Quantum wavefunction engineering of dopant-based Si nanostructures reveals new physics in the 
solid-state, and is expected to play a vital role in future nanoelectronics. Central to any fundamental 
understanding or application is the ability to accurately characterize the deformation of the electron 
wavefunctions in these atom-based structures through electric and magnetic field control. We present 
a method for mapping the subtle changes that occur in the electron wavefunction through the 
measurement of the hyperfine tensor probed by '^^Si impurities. We calculate Stark parameters 
for six shells around the donor. Our results show that detecting the donor electron wavefunction 
deformation is possible with resolution at the sub-Bohr radius level. 

PACS numbers: 71.55.Cn, 03.67.Lx, 71.70.Ej, 85.35.Gv 



The exponential miniaturization of semiconductor 
technology over the past 50 years has ushered in an era 
of nanoscale quantum electronics. At near atomic di- 
mensions, conventional device operations are strongly af- 
fected by quantum phenomena in the solid-state [l|. To 
ensure continued progress in semiconductor electronics, 
and indeed in the drive for new quantum nanoelectronic 
devices, the inherently quantum aspects of such systems 
need to be understood and even incorporated into de- 
vice functionality. The possibility of harnessing quantum 
phenomena in devices has produced revolutionary ways 
of performing computing, as exemplified by the rapidly 
developing fields of quantum computing and spintron- 
ics 0. One central concept of quantum electronics is 
the ability to induce controlled deformations of a specific 
donor-bound electron wavefunction by external electric 
and magnetic fields. Accessing the details of such wave- 
function engineering is critical to understanding and de- 
veloping new devices and applications. However, until 
now there has been no way of quantifying the nature of 
such wavefunction distortions beyond indirect means [l[ . 

In this letter, we propose an electron- nuclear double 
resonance (ENDOR) experiment to directly measure the 
gate induced Stark shift of the donor electron hyperfine 
tensor at specific lattice sites near the donor site (Fig- 
ure 1). ^^Si atoms distributed randomly in the lattice 
provide a direct nuclear spin probe of the donor elec- 
tron wavefunction within the Bohr orbit region. Our 
atomistic tight-binding simulations for lattice regions in- 
volving over a million atoms show that this technique 
provides a spatial map of the bound donor electron re- 
sponse to a controlling gate field to sub-Bohr orbit resolu- 
tion, with excellent correlation to the deformed electronic 
wavefunction, and confirm the feasibility of detecting 
such field induced hyperfine resonance shifts. The tech- 
nique also has wide applicability as it can in principle be 
extended to map out the electric field response of wave- 
functions in single electron Si quantum dots, quantum 
wells or other nanostructures. The ability to map single 



electron wavefunction distortions may have far reaching 
consequences for many current and future quantum na- 
noelectronic applications. 
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FIG. 1: Schematic of the technique. Top row: A series of 
donors in Si under a gate. Inset shows classification of the 
sub-Bohr radii region into shells based on symmetry and dis- 
tance. Bottom row: Probing the field-induced distortions of 
the donor wavefunction by a ^^Si atom using hyperfine inter- 
action. 



Silicon-based quantum nanoelectronic systems bene- 
fit from long spin coherence times and the expertise of 
the semiconductor industry in scalable system design and 
manufacture. As a result there are a number of key pro- 
posals for quantum computing devices, including substi- 
tutional donors 0, 0, [Ejl, 0] , gate-confined 2DEGs 01 , 
and Si quantum well s [9| . Advances in single atom [lO| 
and ion implantation [llj technologies may enable repeat- 
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able fabrication of dopant-based nanostructures. Some 
recent structures include a single gated donor in a Fin- 
FET [H , a gated two donor charge qubit , a 2D gated 



donor layer [l3|, and an ID metallic wire of donors [l4|. 
The wavefunctions of such donor based nanostructures 
vary considerably from their bulk counterparts, yet are 
critical to device operation. A direct map of the wave- 
functions and their electric field response will be impor- 
tant in novel quantum device design and engineering. 

The method we describe here uses the hyperfine inter- 
action between a donor bound electron spin and a nuclear 
spin of a ^^Si isotope in a lattice of spinless ^*Si atoms, 
similar to the method used by Hale and Mieher 15|, but 
critically we include and analyse the effect of a control- 
ling field deforming the donor wavefunction. Although 
current technology limits the substitution of a ^^Si atom 
by a ^^Si atom at a specific point in the lattice, it is 
nevertheless possible to prepare ensemble device samples 
(Fig. 1) with ^^Si atoms distributed randomly around a 
gated donor. The hyperfine interaction between a donor 
electron spin 5* and a ^^Si nuclear spin I is H = I ■ A - S. 
Taking the origin at the ^^Si nucleus, the A tensor is 



f |vl/(0)p-|-(*|5^ 



I*) , (1) 



where 77 and 75 arc the nuclear and electronic gyromag- 



netic ratio respectively, and 



{x, y, z). The first 



term in (1) is the Fermi contact hyperfine interaction, 
denoted here as /3, and is directly proportional to the 
electronic probability density at the ^^Si site. The second 
term represents the magnetic dipolar or anisotropic hy- 
perfine interaction between the two spins, denoted as Bij. 
This dipolar term can also contribute to the ENDOR res- 
onance energies providing a further measure of the dis- 
tribution of the donor electron wave function about the 
probe site. 

ENDOR measurements were performed, first by Feher 
[lit ^iid later by Hale and Mieher [l^ , to study parts of 
the ground state wavefunction of a donor close to the nu- 
cleus, resolving as many as 20 shells (Fig 1) [l^. Several 
theoretical models have calculated the hyperfine tensors 
of a few shells with semi-quantitative agreement with the 
experiments [13, • A recent ab-initio DFT study was 
able to calculate very accurately the tensor components 
of a few shells in the vicinity of the donor nucleus [lit . 
Changes in the Fermi contact hyperfine constants un- 
der a uni-directionally applied stress were also measured 
[iot . It was shown 21 1 that inclusion of the anisotropic 
hyperfine interaction in spin coherence time calculations 
provides remarkable agreement between theory and re- 
cent measurements. The only work on the Stark shift of 
the hyperfine tensors to date computed the Fermi contact 
couphng for 3 lattice sites [il ]. 

The single donor wavefunctions subjected to constant 
electric fields were computed using an atomistic semi- 



empirical tight-binding (TB) model involving a 20 or- 
bital per atom basis of sp'^dPs* (spin) orbitals with near- 
est neigbour interactions. The donor was modeled by a 
Coulomb potential with a cut-off value Uq at the donor 
site. The magnitude of Uq was adjusted to obtain the ex- 
perimental binding energy of Si:P. The total Hamiltonian 
was diagonalized by a parallel Lanczos eigensolver to ob- 
tain the low lying impurity wavefunctions, which were 
then used to evaluate the hyperfine tensors using equa- 
tion (1). The TB method used here is embeded in the 
Nanoelectronic Modeling Tool (NEM0-3D) and has 
been successfully applied to compute Stark shift of the 
donor contact hyperfine coupling (2^ in excellent agree- 
ment with experiments [2^ and with momentum space 
methods [2^. The method was also used to interpret 
single donor transport experiments in FinFETs [l[. 

Two major issues need to be addressed for practical 
implementation of the experiment. First, it is necessary 
to associate each hyperfine resonance peak with a ^^Si 
lattice site. Second, electrostatic gates may give rise to 
inhomogeneous electric fields in the lattice, subjecting 
each donor to a different E-field, and may limit the dis- 
tinguishability of signals corresponding to a ^^Si site. 

The first issue has been resolved in previous works by 
classifying the lattice sites according to their symmetry 
and distance from the donor. For example, the lattice 
sites in the [100] equivalent directions from the donor are 
grouped into a different symmetry class as opposed to 
the points in the [110] or [111] directions. At E = 0, all 
the points equidistant from the donor are responsible for 
a single hyperfine peak, and can be grouped as a shell 
for ease of identification. At a non-zero E-field, these 
points are no longer equivalent, giving rise to multiple 
resonance peaks. As an example, there are 6 points one 
lattice constant qq away from the donor along the [100] 
equivalent directions. If an E-field directed along [010] is 
applied, the four points (with lattice sites (iag, 0, 0) and 
(0, 0, ±ao)) lying in a plane perpendicular to the field are 
still equivalent, and produce a single resonance resonance 
peak. Since the site (0,ao,0) is at a different potential 
than (0, — oo, 0), two separate peaks are observed. Over- 
all, the six-fold degenerate hyperfine peak splits into 3 
components (Fig 2a). 

To minimize the effect of inhomogeneous E-ficlds, it 
is usually advantageous to introduce the donors by ion 
implantation rather than by bulk doping. Ion implan- 
tation at several 100 keVs with typical doses of 10^^ 
cm~^ can ensure a rapidly decaying Gaussian depth dis- 



tribution with a sharply peaked mean depth |25l . 1271 ] . In 
fact, the problem of inhomogeneous E-fields were dealt 
with in this manner in Refs [25[ , which successfully mea- 
sured the Stark shifted hyperfine interaction between the 
donor electronic and nuclear spins. Simulations of the 
ion implantation process can also yield an estimate of 
the typical uncertainty in donor depths, which can be 
incorporated in the data analysis. In addition, clever 
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FIG. 2: Relative change in the contact hyperfine coupling f) 
as a function of electric field for two groups (shells) of points 
along (100) (a, b), {110} (c, d), and (111) (e, f). The coordi- 
nates of the lattice sites are in units of ao/4, where ao = 0.543 
nm. Stark shift of a dipolar component Bzz is also shown in 
(g) and (h) for two shells. The error bars are due to an un- 
certainty of 0.1 MV/m in the E-field. 



gate designs can be utilized that combine the advantages 
of fairly uniform E-fields generated from a parallel plate 
capacitor-like structure, while making sure that the mi- 
crowave radiation from the spin excitations are not com- 
pletely shielded. We have also incorporated reasonable 
uncertainty in the E-ficld in our analysis to predict suit- 
able field regimes for measurements. All these measures 
can help mitigate the effect of inhomogeneous fields. 

In Fig 2, we show the contact hyperfine frequencies as 
a function of E-field for 6 different groups (shells) of sites 
around the donor. Fig 2a and 2b are for sites along [100] 
with distances of ap and 2ao from the donor respectively. 
The degenerate point &i E — Q splits into three curves 
for these groups, as discussed before. Fig 2c and 2d are 
for two shells along [110], while 2e and 2f are for those 
along [111]. The 12 equivalent points in a shell of [110] 
split into 3 groups, while the 4 equivalent points of a shell 
of [111] split into two groups. 

As shown in Fig 2, the frequency axes of the various 
plots do not overlap with each other even at fields of 3 
MV/m. This means that the shells can be distinguished 



TABLE I; Qudratic (772) and linear (771) Stark coefficients for 
the tensor components of 6 shells around the donor in units 
of 10"^ mVMV^ for 772 and 10"^ m/MV for 771. 



13 B,, 
Site 772 771 772 7/1 



Site 



772 7/1 772 771 



(0 4) -3.8 1.7 -3.9 1.8 (0 8) -4.1 1.8 

(0 4 0) -3.5 -6.4 -3.8 4.5 (0 8 0) -2.2 18.7 

(0 4 0) 
(0 2 2) 

(0 2 2) -3.9 5.0 



■3.6 9.8 
■3.8 -1.5 



-3.9 
-3.7 
-3.7 
-3.6 



1.8 
5.8 
-2.4 
-9.7 



-3.7 -1.0 (0 8 0) -2.2 21.5 

-4.3 -1.8 (0 4 4) -3.5 -5.7 

-4.4 5.6 (0 4 4) -3.6 9.1 

(2 2) -3.6 2.2 -4.4 2.0 (4 4) -3.9 3.8 

(1 1 1) -3.8 2.8 -0.1 0.4 (3 3 3) -3.6 -5.6 -0.04 -0.7 

(1 1 1) -3.8 0.7 -0.1 -0.3 (3 3 3) -3.7 9.1 



3.8 13.1 
4.7 3.2 



-0.06 0.8 



even in the presence of electric fields. To show the de- 
tectability of the points within a shell, wc have incorpo- 
rated an uncertainty of 0.1 MV/m (estimated from Fig 
1 of Ref [2^) in the E-field, and represented by error 
bars in Fig 2. This shows that even with some inho- 
mogcneity in the field, most sites can still be identified. 
Furthermore, the distinguishability improves at higher E- 
fields, and for shells farther away from the donor within 
the Bohr-radii due to their increased splitting. The shell 
comprising of the nearest neighbor sites of the donor (Fig 
2e) may not be distinguishable due its small splitting and 
low resonance frequencies. A measurement of the shift 
of the contact hyperfine frequencies directly provides a 
measure of the shift in the electron probability density 
as /?(£;) -/3(0)cx|V'(£;)|2-|V^(0)|2. 

In Fig 2g and 2h, we show the Stark shift of one of 
the dipolar components B^z- These terms are consider- 
ably small in magnitude relative to /3. However, these 
give a measure of the field induced deformation in the 
wavefunction about the donor nucleus through the dipole 
operator (eq (1)). Theoretical models involving donor 
potential and crystal band-structure may be optimized 
to fit these hyperfine components, and thus to improve 
their accuracy. 

We provide quantitative Stark shift data in Table 1 for 
the components (3 and B^z of the six shells considered 
in this work. To concisely present the Stark shifted hy- 
perfine frequencies, we fitted the curves of Fig 2 to the 
form Aa{E) = a{0){r]2E^ + rjiE), where a = {P,B,j), 
rj2 and 771 are the quadratic and the linear Stark coef- 
ficients respectively. The values of 772 and 771 obtained 
from the fit are listed in Table 1. Given the hyperfine 
tensor component aX E = Q for a shell and an applied 
field value, one can calculate both the Fermi contact hy- 
perfine coupling and a dipolar tensor component using 
this table. As an example, the zero-field hyperfine fre- 
quencies of /9 and Bzz, reported in Table II of Ref [13] 
for the first shell along [100], are 2981 and 41.4 kHz re- 
spectively. With E — A MV/m, and using Table I and 
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the qudratic eq, we predict that f3 and B^z of the site 
(0,4,0) decrease by 243 and 3.3 kHz respectively, a net 
change which should be experimentally detectable. In 
comparison, the other two non-equivalent sites, (0,4,0) 
and (0, 0, 4), are shifted by 55 and 160 kHz in /3 and 0.7 
and 2.2 kHz in B^z respectively. Therefore, these 3 sites 
should be distinguishable. 
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FIG. 3: Top row: The stark shifted spectrum of a P donor 
at 3.8 nm from an oxide interface. Bottom row: (a), (b) and 
(c) sliow the P donor ground state wavefunction at tliree dif- 
ferent E-fields, while (d), (e) and (f) show the corresponding 
hyperfine maps in the form of Byy tensor component. The 
E-field is perpendicular to the interface. 

Lastly, we give an illustration of how the method is 
useful for understanding the extent of wavefunction de- 
formation and quantum confinement of direct relevance 
to quantum nanoelectronics Fig. 3 shows the spec- 
trum and the wavefunctions of a P donor at 3.8 nm 
depth from the oxide interface subjected to electric fields 
0, 20 and 40 MV/m (Fig 3a, 3b, and 3c respectively). 
In this regime, the donor wavefunction can be modified 
adiabatically by the field 0, H^, [H, [2^ , as the electron 
makes a transition from a purely Coulomb confined state 
at = to a purely 2D confined state at the inter- 
face at i? = 40 MV/m. In the intermediate field regime 
{E = 20 MV/m), the electron resides in a superposition 
of Coulomb bound and surface bound states (Fig 3b). 
This serves an example of controlled wavefunction engi- 
neering by electric fields. An associated dipolar tensor 
component, Byy for example, is shown on the 2nd row 
of Fig 3 (d, e and f), and refiects the gradual symmetry 
change of the donor wavefunction. 

In conclusion, we proposed the measurement of hyper- 
fine maps of donors as a means of experimentally charac- 
terizing field induced distortions and symmetry changes 



of electron wavefunction. The nuclear spin of a ^^Si atom 
can act as a probe of the donor wavefunction, providing 
a site by site map of electron localization. Such maps can 
help us investigate the unknown electronic wavefunctions 
in novel Si nanostructures for a host of quantum nano- 
electronic applications, and fine tune various modeling 
techniques at the atomic scale. The predictions of the 
Stark shift of the hyperfine tensors for six different shells 
near a P donor indicate that experimental detection of 
engineered wavefunctions is feasible for lattice sites in the 
immediate vicinity of the donor, thus providing a probe 
of the wavefunction at sub-Bohr radius resolution. 
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